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Missions to the outer planets for spacecraft with a solar-electric propulsion system 
(SEPS) and that utilize a single Venus gravity assist are investigated. The trajectories 
maximize the delivered mass to the target planet for a range of flight times. A comparison of 
the trajectory characteristics (delivered mass, launch energy and onboard propulsive 
energy) is made for various Venus gravity assist opportunities. Methods to estimate the 
delivered mass to the outer planets are developed. 
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Nomenclature 


1 c 

= 

linear least square constant 

c 

= 

magnitude of exhaust velocity 

! c 5 

= 

square of the earth-centered hyperbolic excess velocity 

AE 

= 

total propulsive specific energy increment 

AV seps 

= 

velocity increment provided by onboard SEPS 

g 

= 

gravitational acceleration at the surface of the Earth 

Isp * 

= 

specific impulse 

1 

= 

average specific impulse of a SEPS trajectory 

V 

= 

median of I sp * range in all trajectories of this research 

k 

= 

linear least square scale factor 1 

mo, m F 

= 

initial injected mass and delivered mass 

m 

= 

mass flow rate 

£ TOF 


time of flight 


I. Introduction 

I N exploring the outer solar system planets, a planetary gravity-assist has commonly been used since one or more 
gravity-assists have the potential to save propellant, reduce time of flight (TOF), or both. Because of these 
advantages, many previous interplanetary missions (for example, Mariner 10, Voyager I, II, Galileo, Cassini and 
NEAR) exploited the gravity- assist. 1 However, for a ballistic interplanetary mission, there are very little chances that 
all flyby planets are in the desired positions at the desired time. Commonly, the launch may have to be postponed by 
a couple of years if a few-day launch opportunity is missed. Moreover, this short launch opportunity shrinks rapidly 
as the number of flyby planets increases. Solar Electric Propulsion Systems (SEPS) can provide longer launch 
opportunities than a ballistic mission. The trajectory shaping capability of SEPS as well as the longer launch 
opportunity increases mission design flexibility. 
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The flexibility of SEPS missions also allows the gravity-assist to occur at different times even with the same 
TOF and the same launch year constraints. In a previous study, 2 we found that the different gravity- assist timing 
changes the mission performance in terms of delivered mass to target planets in outer-planet, single Venus gravity- 
assist missions. Revolution ratio (R-ratio) was introduced to represent the different Venus gravity-assist timing and 
was defined as the number of Venus revolutions for one revolution of a spacecraft around the Sun. For instance, a 
3:1 R-ratio is one where roughly three Venus years occur during the period from launch to flyby of Venus by the 
spacecraft. For a SEPS mission, the R-ratio can be considered as the problem of trade-off between the launch energy 
and the onboard propulsive energy or similarly, the problem of allocating the flight time to before and after of the 
gravity-assist For this study, we used SEPTOP (Solar Electric Propulsion Trajectory Optimization Program) 3 to 
generate SEPS, outer-planet optimal trajectories. In spite of the necessity of the R-ratio analysis, it can be difficult to 
obtain the different R-ratio optimal trajectories with SEPTOP. In this paper, the R-ratio problem is revisited to 
clarify the underlying fundamentals of the problem and methods to estimate the delivered mass of the different R- 
ratio trajectories are developed to support the R-ratio analysis. 


II. Revolution Ratio 


Figure 1 illustrates the trajectories in an Earth-Venus-Satum (EVS) mission with three different R-ratios. The 
TOF of all trajectories in the figure is 6 years. It is clear that a spacecraft on the trajectory with the largest R-ratio 
spends the most time thrusting before the flyby occurs. Given the nature of SEPS performance - namely, that orbital 
energy addition is more efficient near the Sim due to greater power availability and control authority - it would seem 
that spending more time in close proximity to the Sun before heading outbound toward the destined target would be 
beneficial in delivering more mass. However, a larger R-ratio trajectory typically needs more launch energy, which 
in turn means a larger proportion of the total required HTTP3900, evs, 4 Engine Case 
energy being provided by an inefficient launch 
vehicle rather than the more efficient low-thrust 
engine. Trade-offs between the launch energy and 
the onboard propulsion result in an optimal R-ratio 
of 2:1 for this TOF 6 year mission. Here, the launch 
energy is provided by Delta IV M+(4,2) launch 
vehicle and High Thrust To Power 3900 s (HTTP 
3900) thruster 4 is used as onboard propulsion system. 

Figure 2 shows the delivered mass comparison 
for the EVS, Earth-Venus-Uranus (EVU), Earth- 
Venus-Neptune (EVN), and Earth- Venus-Pluto 
(EVP) missions. This comparison is a result of the 
launch energy and onboard propulsion trade-off. A 
more detailed comparison can be found in Ref. 2. 

Figure 1. Trajectories with three different R-ratios for an 
EVS mission, Delta-IV M+ (4,2). 


6 year TOF ^ — 

Power into thrusters = 23-3825 kW * 
Total Mass Delivered (2:1) = 2493.01 kg 
Total Mass Delivered,<3:l) = 2467.69 kg 
Total Mass Delivered (4:1) = 2226.53 kg 
Propellant Mass (2:1) = 978.69 kg 



III. SEPS Thruster Model 

The free and continuous thrust profile of SEPS 
provides serious difficulties in analytical trajectory 
design. In this paper, a thruster model that represents 
the overall thrust profile is developed to be exploited 
in the delivered mass estimation of SEPS trajectories. 
The thruster model includes constant thrust, mass 
flow rate and or exhaust velocity c . These 
parameters are varying during a mission and the 
variation is mission specific. Because of the variable 
characteristics of a SEPS, the established rocket 
equation is not applicable for a SEPS mission. With 
constant parameters of a SEPS thruster, it would be 
possible to use the rocket equation to estimate the 
delivered mass of the mission. The SEPS thruster 
model is a function of the input power so the 
modeling of the constant thruster parameters is 



Figure 2. Delivered mass of three different R-ratios: EVS, 
EVU, EVN, and EVP missions. 
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equivalent to the selection process of the effective input power (P cf y) within the operating power range. Once the Pcff 
is determined, the thrust and mass flow rate are calculated from the original thruster model. 4 In case the duty cycle 
of the SEPS is not 100%, then the calculated thrust and mass flow rate are reduced by multiplying by the duty cycle 
percentage. For multiple thrusters, the number of thrusters is multiplied to the thrust and mass flow rate to find the 
total thrust and mass flow rate for a mission. With these results, the mission-specific overall 1^ (or magnitude of 
exhaust velocity c) is calculated by Eq. (1). The m is mass flow rate and g is gravitational acceleration at the surface 
of the Earth. 


/ -thrust/ - c/ (1) 

* /mg /g 

If we define 1^* as an average mission-specific I sp , the I sp # is calculated by the initial mass that is inserted to orbit 
by a launch vehicle m 0 , the delivered mass m F and the velocity increment provided by SEPS during the flight AV S eps 
from SEPTOP by Eq. (2) (Ref. 5). 


r * — SEPS (2) 

* g-ln(m F /m 0 ) 


The Isp represents the averaged effect of the thrusting in a trajectory. For all outer-planet trajectories of this 
research, the 1^* varies between 3166.357 s to 3628.537 s. The median of the 1*/ variation I^ 1 " is 3397.4 s and it is 
90% of the maximum 1^ of the HTTP 3900 thruster. The P eff that generate 1^“ of 3397.4 s is calculated from the 
HTTP 3900 thruster characteristics and it results P eff of 4.32 kWe. It is approximately 75% of the maximum 
operating power of the thruster. The I^ 111 (or c) is the representative value for a range of actual I sp ’s so there exists 
differences between the 1^™ and the actual I^’s. The difference can be as large as 230 sec in this research so the 
validity of this modeling strategy should be reviewed before using the I S p m . The validity of the 1^“ as a 
representative value for a SEPS thruster can be shown by examining the sensitivity of the delivered mass (m F ) to the 
variation of I sp . In order to see this, consider the solution to the rocket equation in Eq. (3) (Ref. 5). 



(3) 


Examining first-order variations in I sp and m F . 
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which can be expressed using Eq. (2) as 


•d Vs eps _ dm F 

glsp K m F 

In Eq. (5), AV S EPs^gIsp is typically a number less than one. For a typical EVS mission, AVs^/gl^ is about 3x1 0' 1 . 
For sizable delivered masses of m F = 2000 kg, this makes dm/m F relatively insensitive to dls/I^ in a SEPS mission. 
On the other hand, for a chemical engine mission with the similar AV, the dm p/m F is more sensitive to disp than for a 
SEPS mission because of the smaller I sp range of 200 to 410 s (Ref. 6). 

If the change of 1^ is so large that the first-order analysis is not accurate, one must compare the nonlinear values 
directly. Starting from Eq. (2), let 
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and define 


A m F = m F2 -m F 

= 7 jp2 ~~ 7 apl 
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then 
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and finally. 
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The representative 1^” value for the HTTP 3900 thruster is 3397.447 s and the variation of I sp * in all trajectories with 
various TOF is I^ 111 ± 230 s. The maximum Amj/m F1 deviation for the range of I sp m is less than 5 % for AV S eps = 20 
km/s case. In Fig. 3 the variation of m F is shown for the variation of 1^* near the representative 1^“ value. In the 
figure, the maximum m r estimation error caused by the I sp m modeling in 
acceptable error for preliminary mission 
planning. In summary, the delivered mass is 
not highly sensitive to the specific impulse 
modeling inaccuracy for the SEPS outer- 
planet missions and therefore, a constant 
specific impulse can be used for several 
mission planning purposes. 
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IV, Estimation of Delivered Mass of 
Different R-ratio Trajectories 
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A. Estimation of Delivered Mass with 
Linear Least Square Method 

The bases of the delivered mass 
estimation method are the regularities of the 
relation between the delivered mass and the 
total specific propulsive energy increment. 

The total propulsive specific energy increment is defined as 



SEPTOP ^ W" 

1^=3374.557 s, mps2505 kg 

=3397.447 s, 1^2510.604 kg 
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Figure 3. Delivered mass sensitivity to the specific impulse 
variation. 
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where AE is the total propulsive specific energy increment (excluding an energy increment by a gravity assist) and 
the launch energy C 3 is the square of the earth-centered hyperbolic excess velocity of a spacecraft when it has just 
separated from the launch vehicle. Figure 4 shows the AE of three R-ratio trajectories for EVS to EVP missions. In 


Fig. 4, the AE curves for a target planet have 
similar shape regardless of the different R-ratio. 
In other words, AE curves with different R- 
ratios are related by a shift along with the TOF 
axis and the AE axis. The amount of the 
horizontal shifts (along with the TOF axis) is 
investigated to determine one Venus year 
(224.701 days) and the amount of the vertical 
shifts (along with the AE axis) is found to be 30 
km 2 /s 2 between 3:1 and 2:1 R-ratios and 15 
km 2 /s 2 between 4:1 and 3:1 R-ratios in Fig. 4. 
Figure 5 compares the 2:1 and 4:1 R-ratio AE 
curves from SEFTOP with the shifted 3:1 R- 
ratio AE curve for EVS mission. Overall, the 
shifted curve matches 2:1 and 4:1 R-ratio AE 
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Figure 4. Total propulsive energy increment versus Time of 
flight 



curves accurately. This result means that the AE 
curves with different R-ratio to a target planet 
have similar profiles in spite of their different 
delivered mass profile. There are many 
potential applications of this result; for example, 
it is possible to predict the delivered mass of a 
mission with different R-ratio. Also, the search 
space for optimization can be dramatically 
reduced by previewing the delivered mass 
profiles of the different R-ratio trajectories and 
selecting the most promising R-ratio before 
starting the optimization. The comparison 
results for EVU, EVN and EVP missions are 
not included in this paper for conciseness but 
their results are similar with Fig. 5. 

Another result that should be mentioned is 
that the horizontal distance between the curves 
is constantly one Venus year in all missions. 
Together with the definition of R-ratio, it is an 
interesting result that explains the physical 
differences between the different R-ratio 
trajectories. The difference in TOF from the 
launch date to the flyby date between the 
different R-ratio trajectories is roughly one 
Venus year by the definition of R-ratio so, in 
the point of AE , the effect of the different R- 
ratio is that a very similar AE curve is repeated 
with the TOF difference of one Venus year. 
This result provides a method to generate a 
delivered mass profile of a different R-ratio 
trajectory without actually calculating it by 
SEPTOP reducing time to provide preliminary 
mission planning results. In generating a 
delivered mass profile from a AE curve, the 
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Figure 5. Total propulsive energy increment versus Time of 
flight 
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Figure 6. Delivered mass and AE for 3:1 R-ratio, EVS mission. 


relation between the AE curve and the delivered mass profile needs to be established. 
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Figure 6 shows the delivered mass and AE for 3:1 R-ratio, EVS mission. The delivered mass generation method 
is based on the symmetry between the AE and the delivered mass relative to a horizontal line in Fig. 6. Because of 


Table 1. k and C for EVS to EVP missions 


Target 

R-ratio 

k (kg/km 2 /s 2 ) 

C(kg) 

Saturn 

2:1 

8.501 

3374.258 


3:1 

5.154 

2928.289 


4:1 

3.579 

2481.234 

Uranus 

2:1 

4.908 

2800.503 


3:1 

3.866 

2616.894 


4:1 

2.898 

2329.991 

Neptune 

2:1 

4.163 

2592.702 


3:1 

3.594 

2506.627 


4:1 

2.592 

2269.293 

Pluto 

2:1 

3.964 

2584.694 


3:1 

3.414 

2492.664 


4:1 

2.511 

2233.285 


their symmetry, the delivered mass profile can be represented by AE with a scale factor and a constant. Let//# be 
the delivered mass function and // # be the AE ( # function where £is the time of flight. Then, 

fi(& + k = C (12) 

where k is a scale factor between the delivered mass and AE, and C is the representing constant for the trajectories to 
a target planet with the R-ratio. Equation (12) is rearranged to be used in linear least square fitting in the below 
equation. 

k-M&-c=-f,(0 (13) 

For the given TOF * = A-**, n, Eq. (13) can be used to construct below linear least square equation. 


7(6) -i' 
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The Q-R factorization 6 7 is used in solving Eq 
to EVP missions, the k's and Cs are calculated 
and shown in Table 1. 

The thruster used is HTTP 3900 and the launch 
vehicle is Delta IV M+(4,2). The k and C 
depend on the thruster and the launch vehicle 
because the relation between the AE and the 
delivered mass is determined by them. The k 
represents the propulsion efficiency of AE to 
the delivered mass. Large k means more 
efficient thrusting (launch vehicle + SEPS) for 
a given delivered mass. The unit of k is kg ■ s 2 / 
km 2 . The constant C can be interpreted as a 
limit of the maximum delivered mass of the R- 
ratio trajectory class to a target planet. 
According to Table 1, k ' s are larger for closer 
targets and that can be interpreted to mean that 
the thrusting is more efficient than for a farther 
targets. At the same time, Cs for farther targets 
are smaller than Cs for closer targets so the 


. (14) to calculate the k and C. From AE ' s of three R-ratios for EVS 



Figure 7. Delivered mass comparison between SEPTOP and k> 
C fitting for EVS mission. 
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limit of the delivered masses for farther targets 
are smaller. In Eq. (12), the k and C are 
regulating each other in the sense of increasing 
or decreasing the delivered mass to a target 
planet. The delivered mass is determined by the 
interplay between the k and C and the given 
profile of AE as a function of TOF. 

In order to check the validity of k and C, the 
delivered mass profiles are regenerated using 
Eq. (13) and k and C of Table 1. Figure 7 shows 
the comparison of the delivered masses from 
SEPTOP with the delivered mass from K C and 
AE from SEPTOP for an EVS mission. The 
overall predictive performance is observed in 
the figure. These results mean that if the k , C 
and AE were known for a mission, the delivered 
mass profile for the mission could be estimated 
without using SEPTOP, As previously 
explained, the AE curves for a target planet can 
be approximated by a planar shift of a known 
AE curve. At present, a AE curve to a target 
planet is required to be calculated by SEPTOP 
to approximate the AE of the different R-ratio 
trajectories to the target planet. On the other 
hand, the k and C for a mission can be closely 
approximated once the k and C to a different 
target planet but with the same R-ratio are 
known from the SEPTOP result. In other words, 
k and C of a mission to a target planet with an 
R-ratio have a relation with the k and C of a 
mission to a different target planet with the 
same R-ratio. The relation between it’s and Cs 
with the different R-ratio to a target planet 
needs more research to be clarified. The 
relations between fc’s and Cs of different 
missions are shown in Figs. 8 and 9 as 
functions of inverse of the semimajor axis of 
target planets. In the figures, the k's and Cs are 
shown as linear functions of the inverse of 
semimajor axis. It is also possible to calculate 
&’s and Cs for other target planets if one use 
the linear functions and then possible to 
estimate the delivered mass to the target planets. 
In summary. Fig. 10 illustrates a procedure of 
linear least square method. 

Figure 11 and Fig. 12 shows the 
comparisons between the estimated delivered 
mass and the delivered mass from SEPTOP for 



Figure 8. k variation as a function of inverse of semimajor axes. 



Figure 9. C variation as a function of inverse of semimajor axes. 



Figure 10. 
method. 


A procedure of linear least square estimation 


EVS and EVP missions. The estimated delivered mass is calculated with k and C and planar shifted AE curves. EVU 
and EVN missions show similar results so they are excluded for the conciseness of the paper. 

In Figs. 1 1 and 12, there are some differences between the regenerated delivered mass and the SEPTOP result in 
spite of the overall coincidences. The differences in short TOF range are caused by the incomplete symmetry 
between the delivered mass profile and the AE profile. The incompleteness is more outstanding in all of the 2: 1 R- 
ratio cases than the other cases and this is the reason of the relatively larger differences in short TOF, 2:1 R-ratio 
cases. The source of the differences in intermediate to long TOF range is the linear least square error. The linear 
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least square algorithm tries to minimize the 
total error in all of the data range so a small 
range of large error usually degrade the total 
accuracy and the error spread out over all the 
data range. Therefore, there are relatively 
larger errors in long TOF, 2:1 R-ratio cases 
than in other R-ratio cases. The maximum 
difference between the regenerated delivered 
mass and the SEPTOP result is about 500 kg 
in short TOF. This estimated delivered mass 
can be used for some mission planning 
purposes, because the estimated delivered 
mass curve is a good approximation for 
middle to long TOF range missions, however, 
a more precise method to predict the 
delivered mass is required to determine 
which R-ratio trajectory is the best candidate 
for a given TOF range of interest without 
spending lots of effort in SEPTOP. 



Time of Right (year) 


Figure 11. Delivered mass comparison, the estimated versus 
SEPTOP, EVS mission. 
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B. Estimation of Delivered Mass with 
Hybrid Method 

The linear least square estimation method 
has some amount of error especially in its 
short TOF ranges. The error is caused by the 
incomplete symmetry between the delivered 
mass profile and the AE profile. In spite of 
the incomplete symmetry, however, the AE 
profile has significant advantage because the 
similarity in AE profiles among the different 
R-ratio cases provides a basis to predict the 
delivered masses of the different R-ratio 
cases. A more precise estimation method, a 
hybrid method, also relies on the similarity of 
AE profiles as the linear least square method 
but it also employs a method that uses the 
thruster model and rocket equation. 

In the hybrid method, the TOF range is 
divided to two ranges, short TOF and remaining TOF range, and a method that uses the rocket equation is applied to 
that short TOF range whereas the linear least square method that uses k and C is applied to the remaining TOF range. 
First, a solution to the rocket equation is shown below. 5 



Time of Rigtrt (year) 

Figure 12. Delivered mass comparison, the estimated versus 
SEPTOP, EVP mission. 


m F 


- m 0 exp 



(15) 


The delivered mass can be calculated given m#, AVseps and c. The way to determine the m 0 and AVseps for a given 
TOF is explained below and the c is modeled by a constant l sp SEPS model. 

In estimating the delivered mass, the given data are the AE and the delivered mass profile of an R-ratio case (3:1 
R-ratio case in this research). By a planar shift, the AE profile of 2:1 and 4:1 R-ratio cases can be approximated. 
However, it is required to find the AV S eps out of AE because Eq. (15) needs AVseps to calculate the m F and AE is the 
sum of the C? and the square of AV S eps * Of course, the AVseps and also C 3 are assumed to be known for an R-ratio 
case in order to calculate the AE profile that is to be shifted. Therefore, the C/AE profiles of all R-ratio cases are 
compared to find the relation between them. Figure 13 shows the variation of C/AE versus AE profile for the short 
TOF range, for all target missions. By comparing the AE values with Fig. 4, one can notice that the TOF gets shorter 
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as AE value increases in Fig. 13. The C/AE 
profiles are generated with SEPTOP results 
and almost identical in all R-ratio missions. 

This fact provides a way to find the AV S eps 
out of the shifted AE profile. If we assume 
the C/AE profiles of other R-ratio cases are 
identical with the C/AE profile of the known 
R-ratio case, which is true for this research, 
the Cj and AV SE ps can be extracted from the 
shifted AE profile and the extracted AV S eps 
can be put into the rocket equation solution to 
predict the delivered mass of other R-ratio 

cases. The reason why the C/AE profiles are 
almost identical in all missions can be Figure 13. C 3 /AE vs. AE, short TOF, all R-ratios, all target 

explained as follow. Because all missions use missi0ns * 

the same launch vehicle. Delta IV M+(4,2) and the same SEPS (HTTP 3900 thruster), and the shapes of trajectories 
are similar except its R-ratio, the portion of C 3 (or in AE profile, i.e. C/AE profile, are identical within an 

allowable error range for all target missions. Figure 14 illustrates the procedure of the hybrid estimation method for 
the short TOF range. Once we have a trajectory data for an R-ratio case from SEPTOP, the AE profile and the C/AE 
profile of the R-ratio case are available for the next step. By shifting the AE profile in TOF -AE plane, the AE profile 
of different R-ratio case can be approximated. The AV SEPS is computed with the assumed C/AE profile and the 



definition of AE. Then with the constant 
exhaust velocity from the constant thruster 
model and the AV S eps* the delivered mass is 
computed for a different R-ratio case without 
using SEPTOP. If the trajectories in different 
missions significantly differ from each other, 
then this assumption in the similarity of 
C/AE profile is not valid. In this research, 
however, all trajectories have similar 
characteristics except the different R-ratios in 
order to compare the effect of R-ratio to the 
delivered mass performance of a trajectory. 

Another important application of the 
C/AE profiles is its peak point. The peak 
point is interpreted as the maximum relative 
usage of C 3 and there is a change in thrusting 
strategy in the left and right parts of the peak 
point. The peak point is equivalent to the 
beginning of stationary state of AE profile in 
Fig. 4. 

Because the C/AE profiles of other R- 
ratio cases are assumed to be identical with 
the known C/AE profile, the peak C/AE 
points of other R-ratio cases are also identical 
with the known peak point. The peak point is 
used as the border between the short TOF 
range and the remaining TOF range in the 
hybrid method. Therefore, the method that 
uses the rocket equation solution is applied 
for the TOF < peak TOF (the TOF value for 
the peak point) and the previous method that 
uses k and C is applied for the TOF > peak 
TOF. 

The m 0 in Eq. (15) can be calculated by a 



Figure 14. Procedure of hybrid method for short TOF range 
estimation. 




Figure 15. Hybrid Method Delivered mass comparison, the 
estimated versus SEPTOP, EVS mission. 


9 

American Institute of Aeronautics and Astronautics 










given launch vehicle model if Cj is known. The 
C s and AV S eps can be extracted from the 
assumed C/AE profile and the shifted AE 
profile as previously explained. The exhaust 
velocity c is not a constant for a SEP S but it is 
modeled as a constant with sufficient accuracy. 
Now, one can use the rocket equation solution 
to predict the delivered mass for the short TOF 
missions. The AE of 3:1 R-ratio case is chosen 
to be shifted and its C/AE is used for all other 
R-ratio cases. Figure 15 and Fig. 16 show the 
comparison between the estimated delivered 
mass curves and the delivered mass from 
SEPTOP for all R-ratio cases for EVS and EVP 
missions. The estimated delivered mass curves 
consist of two parts and their border is the peak 
TOF point of C/AE profile. It is clear that there 
is significant enhancement of the accuracy of 
the delivered mass estimation compared with 


3000 


2000 


1000 




Figure 16. Hybrid Method Delivered mass comparison, the 
estimated versus SEPTOP, EVP mission. 


Figs. 10 and 1 1 of the linear least square method. The maximum error in the short TOF range is less than 70 kg and 
the maximum error in the long TOF range is less than 100 kg for the hybrid method. For a comparison, the 
maximum error in the linear least square method was about 500 kg. The reasons of this accuracy enhancement are 
first, the similarity of the C/AE profiles in short TOF range of the different R-ratio cases and second, the better 
linear least square fitting with smaller TOF range for the long TOF range. On the other hand, the error of the hybrid 
method comes from the assumptions about the AE profile, C/AE profile and the modeling error of the constant 
exhaust velocity of SEPS for the short TOF range. The incomplete symmetry between the delivered mass profile and 
the AE profile and the linear least square error are the sources of the errors in the longer TOF range. 


V. Conclusion 

Various R-ratio trajectories to outer planets are investigated in this paper. The trade-off between the launch 
energy and the onboard propulsive energy of SEPS for Venus gravity-assist missions is explained. An estimation 
method is presented for predicting mission performance for different R-ratio trajectories and can successfully 
regenerate the delivered mass profile only from one converged set of mission data and a given thruster model. 
Extensions of this method are currently being investigated. 
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